import ij.ImagePlus;
import ij.plugin.PlugIn;
import ij.process.ColorProcessor;
import ij.process.ImageProcessor;


import java.awt.Color;


import java.awt.Font;

import java.awt.Polygon;
import java.io.PrintStream;

import java.util.Locale;
import java.util.Random; 



import pluginUtils.Circle;
import pluginUtils.CircleArrangement;
import pluginUtils.Complexe; 
import pluginUtils.DoyleCircle;
import pluginUtils.Droite;
import pluginUtils.MobiusTransform;
import pluginUtils.CircleInversion;
import pluginUtils.PackCircle;
import pluginUtils.Point2D;
import pluginUtils.Pos3D;
import pluginUtils.Sphere;
import pluginUtils.SphereInversion;
import pluginUtils.Vertex;


public class Lunule extends ImagePlus implements PlugIn{
	{Locale.setDefault(Locale.US); }
	
	private ImageProcessor primeIP; 
	private int X=2000,Y=2000; // image size
	
	private PrintStream output;
	private double xShift=0; 
	private double zShift=0; 
	private double limit,lmoins; 
	
	

	private Random generator=new Random(); // One always needs some randomness
	
	// For the scaling of the image
	private double x0=0,x1=1,y0=0,y1=1;
	private double acoef,ccoef,vcoef,wcoef; 
	
	private String[] palette={"Yellow","Red","Blue","Green","White","NavyBlue","Turquoise","Orange"};
	private Color[] palette1={Color.YELLOW,Color.RED,Color.BLUE,Color.GREEN,Color.MAGENTA,Color.ORANGE};
	public void setCoefs(){
		acoef=X/(x1-x0); 
		vcoef=Y/(y1-y0); 
		ccoef=-acoef*x0; 
		wcoef=-vcoef*y0; 
	}
	 
	
	
	public void setBounds(double xl,double xu,double yl,double yu){
		x0=xl; 
		x1=xu;
		y0=yl;
		y1=yu; 
	}
	
	private void drawCircle(PackCircle c){
		//if(intersect(c))return; 
		
		//if(c.getRadius()<0.1) return; 
		int intX=(int)(c.getX()*acoef+ccoef);
		int intY=(int)(c.getY()*vcoef+wcoef);
		int intRx=(int)(c.getRadius()*acoef);
		int intRy=(int)(c.getRadius()*vcoef);
		int rank=((int)(c.getCouleur()*60))%6; 
	
			
			double Xs=c.getX(); 
			double Ys=c.getY(); 
			double rapp=Xs*Xs+Ys*Ys+1; 
			double u=Xs/rapp; 
			double v=Ys/rapp; 
			double w=(Xs*Xs+Ys*Ys)/rapp;
			
			double Xsr=c.getX(); 
			double Ysr=c.getY()+c.getRadius(); 
			double rappr=Xsr*Xsr+Ysr*Ysr+1; 
			double ur=Xsr/rappr; 
			double vr=Ysr/rappr; 
			double wr=(Xsr*Xsr+Ysr*Ysr)/rappr;
		
			double nrad=(u-ur)*(u-ur)+(v-vr)*(v-vr)+(w-wr)*(w-wr);
			//nrad=Math.sqrt(nrad)*1.414; 
			
			String s="sphere{<"+u+","+v+","+w+">,";
			
			
			s+=nrad; 
			s+=" texture{Texture"+rank+"} finish{Finish"+rank+"}}"; 
			output.println(s);
		
	}
	private void drawCircle(Point2D C,double r){
		double xCenter=C.getX(); 
		double yCenter=C.getY(); 
	 
		
		int intX=(int)((xCenter)*acoef+ccoef);
		int intY=Y-(int)((yCenter)*vcoef+wcoef);
		int intRx=(int)(r*acoef);
		int intRy=(int)(r*vcoef);
		
		
		primeIP.drawOval(intX-intRx,intY-intRy,intRx*2,intRy*2);
	}
	
	private void drawColoredCircle(PackCircle c){
		double xCenter=c.getX(); 
		double yCenter=c.getY(); 
		//if(c.getRadius()<0.001) return ; 
		
		int intX=(int)((xCenter-xShift)*acoef+ccoef);
		int intY=(int)((yCenter-zShift)*vcoef+wcoef);
		int intRx=(int)(c.getRadius()*acoef);
		int intRy=(int)(c.getRadius()*vcoef);
		
		
		primeIP.drawOval(intX-intRx,intY-intRy,intRx*2,intRy*2);
		
		/*
		primeIP.moveTo(intX,intY); 
		primeIP.drawString(c.getMyOrder()+"");
		*/
		
		String s="sphere{<"+xCenter+",0,"+yCenter+">,";
		s+=c.getRadius(); 
		s+=" scale 0.4*y texture{Texture"+generator.nextInt(7)+"} finish{Finish0}}"; 
		System.out.println(s); 
		output.println(s);
	
	}
	private void drawPoint(double x,double y){
		int intX=(int)(x*acoef+ccoef);
		int intY=Y-(int)(y*vcoef+wcoef);
		primeIP.fillOval(intX-10,intY-10,20,20);
	}
	
	private void drawPoint(Point2D A){
		drawPoint(A.getX(),A.getY()); 
	}
	
	private void drawTriangle(double[] X,double[] Y){
		int xint[]=new int[3]; 
		int yint[]=new int[3];
		for(int i=0;i<3;i++){
			xint[i]=(int)(X[i]*acoef+ccoef); 
			yint[i]=(int)(Y[i]*vcoef+wcoef);
		}
		primeIP.fillPolygon(new Polygon(xint,yint,3)); 
		
	}
	
	private void drawTriangleEqui(double[] X,double[] Y){
		double XP[]=new double[3]; 
		double YP[]=new double[3];
		for(int i=0;i<3;i++){

			XP[i]=(X[i]+Y[i]/2);
			YP[i]=(Math.sqrt(3)/2*Y[i]);

	/*
			XP[i]=(X[i]+Y[i]/2); ///Math.sqrt(2);
			YP[i]=(Math.sqrt(3)/2*Y[i]); ///Math.sqrt(2);
			*/

		}
		drawTriangle(XP,YP); 
		
	}
	
	
	
	private void drawLine(double xd,double yd,double xf,double yf){
		int XD=(int)(xd*acoef+ccoef); 
		int YD=(int)(yd*vcoef+wcoef); 
		int XF=(int)(xf*acoef+ccoef); 
		int YF=(int)(yf*vcoef+wcoef);
		primeIP.drawLine(XD,Y-YD,XF,Y-YF); 
	}
	
	private void drawLine(Point2D A,Point2D B){
		drawLine(A.getX(),A.getY(),B.getX(),B.getY()); 
	}
		

	
	public void drawDroite(Droite d){
		double y1=d.getY(limit);
		double y2=d.getY(lmoins); 
		double x1=d.getX(limit); 
		double x2=d.getX(lmoins);
		Point2D p1[]=new Point2D[2]; 
		int indice=0; 
		if((y1>=lmoins)&&(y1<=limit))
			p1[indice++]=new Point2D(limit,y1); 
		if((y2>=lmoins)&&(y2<=limit))
			p1[indice++]=new Point2D(lmoins,y2); 
		if((x1>=lmoins)&&(x1<=limit))
			p1[indice++]=new Point2D(x1,limit); 
		if((x2>=lmoins)&&(x2<=limit))
			p1[indice++]=new Point2D(x2,lmoins);
		
		drawLine(p1[0],p1[1]); 
		
	}
	
	
	
	public void run(String arg) {
		primeIP = new ColorProcessor(X,Y);
		double r2=Math.sqrt(2);
		
		 limit=2.5; 
		 lmoins=-2.5;
	
		setBounds(lmoins,limit,lmoins,limit); 
		setCoefs(); 
		
		
		
		
		
		
		Polygon background=new Polygon(new int[]{0,0,X,Y},new int[]{0,X,Y,0},4);
		primeIP.setColor(Color.WHITE);
		
		primeIP.fillPolygon(background);
		
		
		primeIP.setFont(new Font(null,Font.ITALIC,60));
		primeIP.setColor(Color.BLACK);
		primeIP.setLineWidth(3); 
		
		drawLine(lmoins,0,limit,0); 
		drawLine(0,lmoins,0,limit); 
		setProcessor("Verifying", primeIP);
		
		int nbDiv=50;
		
		double R=1.2; 
		double l=0.9; 
		
		
		Point2D[] lesPoints=new Point2D[2*nbDiv+2]; 
		lesPoints[0]=new Point2D(0,Math.sqrt(R*R-l*l/4));
		System.out.println(lesPoints[0]);
		lesPoints[2*nbDiv+1]=new Point2D(0,-Math.sqrt(R*R-l*l/4));
		drawPoint(lesPoints[2*nbDiv+1]);
		for(int i=1;i<=nbDiv;i++){
			double ti=Math.tan(Math.PI/2-Math.PI*i/(nbDiv+1.0));
			double delta=Math.sqrt(4*R*R*(1+ti*ti)-ti*ti*l*l); 
			double x=(-l+delta)/(2*(1+ti*ti));
			
			lesPoints[2*i]=new Point2D(x,ti*x);
			x=(l+delta)/(2*(1+ti*ti));
			lesPoints[2*i-1]=new Point2D(x,ti*x);
		} // for
		drawLine(lesPoints[0],lesPoints[1]);
		drawLine(lesPoints[0],lesPoints[2]);
		System.out.println("ok"); 
		drawLine(lesPoints[2*nbDiv+1],lesPoints[2*nbDiv]);
		drawLine(lesPoints[2*nbDiv+1],lesPoints[2*nbDiv-1]);
		for(int i=1;i<nbDiv;i++){
			
			drawLine(lesPoints[2*i],lesPoints[2*i+2]);	
			
			drawLine(lesPoints[2*i-1],lesPoints[2*i+1]);	
		}// for
		
		 show();

	}	
}
	